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Abstract. Attenuation of high-energy gamma rays by pair-production with UV, optical and IR background photons provides 
a hnk between the history of galaxy formation and high-energy astrophysics. We present results from our latest semi-analytic 
models (SAMs), based upon a ACDM hierarchical structural formation scenario and employing all ingredients thought to 
be important to galaxy formation and evolution, as well as reprocessing of starlight by dust to mid- and far-IR wavelengths. 
Our models also use results from recent hydrodynamic galaxy merger simulations. These latest SAMs are successful in 
reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color 
bimodality. We have created 2 models that bracket the likely ranges of galaxy emissivities, and for each of these we show 
how the optical depth from pair-production is affected by redshift and gamma-ray energy. We conclude with a discussion 
of the implications of our work, and how the burgeoning science of gamma-ray astronomy will continue to help constrain 
cosmology. 



1. INTRODUCTION 

The extragalactic background light (EBL) consists of 
photons emitted at all epochs by stars and AGN, and sub- 
sequently modified by redshifting and dilution due to the 
expansion of the universe. The bulk of the EBL energy is 
in two peaks in wavelengths, one in the optical and near- 
IR consisting of direct light from stars, and another in 
the far-IR producing by thermal radiation from dust (see 
Figure 1). At still longer wavelengths, the background 
spectrum is dominated by the cosmic microwave back- 
ground. Emissions from polycyclic aromatic hydrocar- 
bon (PAH) molecules also form an important component 
at a fixed series of wavelengths in the mid-IR. 

Because the production of the EBL is directly linked 
to the star formation history of the universe, limits on the 
EBL can be used to provide constraints on the history 
of galaxy formation and evolution. Measurement of this 
light through direct observation is complicated by fore- 
ground emission from our own galaxy and reflected zo- 
diacal light from our sun, which are much brighter than 
the EBL across most of the optical and IR spectrum [1]. 
Interplanetary dust is the major contributor of foreground 
light at most wavelengths, with starlight becoming sub- 
stantial in the optical and near-IR, and the interstellar 
medium most important in the submillimeter regime. In a 
series of papers by Rebecca Bernstein [2, 3, 4] the optical 
EBL was estimated by measuring the total background 
in 3 HST WFPC2 bands, and subtracting away the zo- 
diacal and galactic contributions. The near-IR flux has 
been measured by the DIRBE instrument on the Cosmic 
Background Explorer [5, 6, 7, 8, 9], using foreground 
source subtraction techniques and modeling of the zo- 
diacal light, and yielding high estimates in this range. 



Results from IRTS [10] claimed an even higher back- 
ground level below 2 jim a level which would require 
a large contribution from an as-of-yet undetected source 
type. An analysis of the fluctuations observed by DIRBE 
[11] yielded upper limits to the EBL in the mid to far- 
IR, but this author's interpretation of a high background 
arising from sources at z > 8 has been challenged by the 
results of other studies [12, 13]. In the far-IR regime, ob- 
servations from the COBE FIRAS [14] and DIRBE [15] 
instruments provide direct measurements of the back- 
ground radiation. At shorter wavelengths (60 and 100 
jim) the DIRBE instrument places upper limits on the 
background. Analyses of the DIRBE data at far-IR wave- 
lengths have been attempted by other authors using new 
models for the zodiacal light and instrument calibra- 
tion [16, 17, 18]. The rather high detection claimed by 
Finkbeiner et al. at 60 and 100 /im, which required care- 
ful subtraction of the bright foreground, has been dis- 
puted by other authors [19]. 

Integration of galaxy counts is a way to set firm lower 
limits on the EBL, although the degree to which these 
measurements converge on the true value remains con- 
troversial. Measured luminosities of galaxies in the Hub- 
ble Deep Field (HDF) [20] place a lower bound on EBL 
fluxes in the 0.4-1 /xm range, a much lower level than 
that proposed in the Bernstein papers, with 2MASS data 
extending counts in this paper to 3 /im. It is argued in 
this work that the flatness of the faint counts in all bands 
(flatter than 0.4 on a log-counts vs. magnitude diagram) 
indicate convergence, and therefore a low EBL across the 
optical and near-IR. A similar analysis conducted using 
IR data from the Subaru deep field [21] modeled possi- 
ble selection effects and supported the conclusion of a 
low optical-near-IR background compared to direct de- 



tection claims. In the UV, limits exist from GALEX [22] 
and observations of the HDF with the STIS instrument 
[23], with the latter finding a considerably higher bound. 
The IRAC instrument on Spitzer has placed lower lim- 
its on several bands in the near- to mid-IR [24], which 
similarly are well below the IRTS and DIRBE direct 
detection fluxes. In a complementary approach to their 
DIRBE sky-subtraction papers, Levenson and Wright 
used IRAC data to calculated the best-fit flux at 3.6 jj.m 
using a profile-fit to estimate the light from the unob- 
servable faint fringes of galaxies, and a broken power 
law model for the number count distribution [25]. ISO- 
CAM and the MIPS instrument on Spitzer have reported 
lower limits from number counts at 15, 24, 70, and 160 
|im [26, 27, 28, 29, 30]. Galaxy counts from the SCUBA 
instrument provide a lower limit at 850 microns [31], 
with an estimated 20 to 30% of the background at these 
wavelengths resolved into point sources. Our previous 
EBL model [32] was systematically low in the mid- and 
far-IR, and we warned in that publication that the model 
should not be trusted in this regime. Our new models, 
presented here, do a much better job in matching con- 
straints from number counts. 

The EBL provides an important test of star-formation 
observations; if there is a large gap between the flux from 
sources in deep surveys and direct measurements, then 
an unresolved population must exist to provide the miss- 
ing energy. The present day EBL is Unked to the history 
of star formation through the stellar initial mass function 
(IMF), as well as the evolving metallicity and dust distri- 
bution, as we will discuss later. It is well established that 
star formation rates per unit cosmological volume were 
much higher in the past, peaking at roughly z^2 [39, 40]. 

Studies of the buildup of stellar mass have typically 
found that the integrated star formation rate tends to ex- 
ceed stellar mass as measured by IR (fossil) light from 
old stars, with the discrepancy becoming worse with red- 
shift. The local star-formation rate is typically measured 
from the strength of the Ha (Balmer) line. At higher red- 
shift, other tracers must be used. These can include UV 
continuum emission corrected for dust, radio emission 
by supemovae remnants, emission lines which are ex- 
cited by UV radiation from young stars, or measuring 
reprocessed light from dust in the infrared. All of these 
methods are fraught with uncertain amounts of contam- 
ination from quasars and biases from dust and metal- 
licity, and different methods can produce widely vary- 
ing results for the same galaxy. An overview of these 
issues can be found in [41]. One solution to the SFR- 
fossil mass discrepancy may lie in a variable IMF which 
suppresses formation of low mass stars in starbursts oc- 
curring primarily at higher redshift [42, 43]. In this way, 
high mass stars that produce tracers of the star-formation 
rate could be produced at a high rate while lower mass 
stars that account for most of the integrated steUar mass 



are created in smaller quantities. Alternatively, it is pos- 
sible that systematic biases in determinations of the star 
formation rate at increasing redshift are responsible for 
the discrepancy. This can be phrased in terms of mea- 
surements of Dave's star-formation activity parameter, 
which determines the Hubble times required for galaxy 
to reach their current mass at current SFR, and is seen to 
decrease too much between present day and z^2. This 
conflicts with model predictions, and the passive popu- 
lation that would be required to counterbalance the high 
star formation in these rapidly growing systems is not 
observed. The recent work of [44] has found lower spe- 
cific star formation rates in galaxies from the SDSS and 
DEEP2 surveys than typicaUy measured by more usual 
means using an alternative method based on higher or- 
der Balmer lines. The lower normalization found in this 
model can be attributed to a number of factors, though 
the possibiUty of a changing IMF cannot be ruled out. 

Unresolved primordial (population III) stars could in 
principle provide a large contribution to the background 
light, and the near-IR peak observed by the IRTS satel- 
lite [10] has been interpreted as the redshifted photons 
from massive population-Ill stars beyond redshift ^9 
[45]. This peak reaches nearly a factor of 10 higher than 
levels from resolved number counts, but this claim is 
at odds with [20, 21, 12, 13] in which it is argued that 
the contribution from resolved galaxies has nearly con- 
verged, and there is little room for any additional sources 
such as an early generation of stars with a top-heavy 
IMF. The levels of metal formation from massive stars 
in these models required metals to either be locked away 
in compact products or dispersed in the IGM in a very 
inhomogeneous fashion. The large fraction of baryons 
which would have to be processed through primordial 
stars, ~10%, and the lack of J-band dropout detections 
of these sources have strongly disfavored this interpre- 
tation [46, 47], as does amount of material that would 
be contained in intermediate mass black holes [48]. Ad- 
ditionally, an extragalactic origin for this feature would 
have a huge impact on the optical depth of TeV gamma 
rays, as we wiU discuss shortly. 

Three approaches have been followed to calculate the 
EBL (as described by [49]): (A) Evolution Inferred from 
Observations - e.g., [Kneiske et al. 2002, Franceschini 
et al. 2008]; (B) Backward Evolution, which starts with 
the existing galaxy population and evolves it backward in 
time - e.g., [50]; and (C) Forward Evolution, which be- 
gins with cosmological initial conditions and models gas 
cooling, formation of galaxies including stars and AGN, 
feedback from these phenomena, and light absorption 
and re-emission by dust - e.g., [51, 32]. We will concen- 
trate on Forward Evolution here, based on semi-analytic 
models (SAMs). 




FIGURE 1. The predicted z=0 EBL spectrum from our fiducial (blaclc) and low (dashed blue) models, compared with experiments 
at a number of wavelengths. Our previous model [32] is also shown for comparison (dotted curve). We have also shown the model 
of Franceschini et al. [33] for comparison (dash-dotted orange line). Upward pointing arrows indicate lower bounds from number 
counts; other symbols are results from direct detection experiments. Some points have been offset slightly to improve readability. 
Lower limits: The open blue-violet triangles are results from Hubble and STIS [23], while the magenta triangles are from GALEX 
[22]. The green and red triangles from Hubble Deep Field [20] and Ultra Deep Field [34] respectively, combined with ground 
based-data. Gold triangles are also from this same work by Madau & Pozzetti. Open red triangles are from IRAC on Spitzer [24], 
and the pink point at 15 flm is ISOCAM [26] on ISO. The remaining lower limits are from MIPS at 24, 70, and 160 /Jm on Spitzer 
[29, 27, 28, 30]. Direct Detection: The higher open blue triangles are from Bernstein [4], while the lower points are the original 
determinations from [2]. The high-reaching cyan points in the near-IR are from [10]. While the rest of the points in this region are 
based upon DIRBE data with foreground subtraction [6] (dark red squares), [8] (orange 4-stars), [25] (red diamond), [7] (purple 
open hexes), [5] (green square), [9] (red asterisks). In the far-IR, direct detection data is shown from DIRBE [18] (blue stars), and 
[15] (green stars), and also purple bars showing the detection of FIRAS [14]. 
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FIGURE 2. The z=0 luminosity density of our model. As 
before, the fiducial model is shown as a solid black line, the 
low model as dashed blue. Data at a number of wavelengths is 
shown from GALEX, SDSS, 6dF [35], 2MASS [36, 37], and 
IRAS [38]. 

2. SEMI-ANALYTIC MODELING 

The use of semi-analytic models to study galaxy forma- 
tion was pioneered by [60], and initial work was carried 
out primarily by two groups based in Durham and Mu- 
nich. Our group was the first to develop a complete SAM 
code for galaxy formation [61] and use it to model the 
EBL. Our code used a merger-tree formalism which was 
consistent with high-resolution dark matter simulations 
of the time, and several models of star formation and su- 
pernovae feedback were investigated. In [62], these mod- 
els were combined with the improved dust emission code 
of [63, 64] for a prediction of the EBL and absorption 
of gamma rays [51]. At the 2004 Symposium on High 
Energy Gamma Ray Astronomy our group presented a 
modified spectrum and attenuation predictions [32]. This 
SED was somewhat lower than our 2001 result [51], pri- 
marily because our 2005 model [32] had been recali- 
brated to fit the local luminosity density as determined 
by new surveys such as 2MASS [36], SDSS [65], and 
2dF [66]. While the prediction of this model was con- 
sistent with optical and near-IR number counts data, it 
was below both the ISOCAM lower limit at 15 microns 
and the direct detection fluxes of DIRBE and FIRAS in 
the far-IR, indicating that the light remission by dust was 
being significantly underrepresented. 

Our current SAM [52] is normalized using a variety 
of updated survey data, particularly in the IR, and also 
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FIGURE 3. The star formation rate density for the fiducial 
(blue) and low (broken orange) models; this plot is from [52] 
(Figure 14). Lower curves of the same colors indicate the 
contribution from starbursts. The data in this plot is from the 
compilation of [39] and has been converted to a Chabrier initial 
mass function. The gray line is the best fit to the data from 
[41]. The thin gray broken contours are estimates from GALEX 
observations [53] 
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FIGURE 4. The integrated stellar mass vs redshift in each of 
our models, also from [52] (Figure 14). As in the previous fig- 
ure, the solid blue curve denotes the fiducial model, and broken 
orange the low model. The solid square is the z=0 estimate of 
[37], the circles are from [54], and the open diamonds are from 
the COMBO-I7 estimates of [55]. The gray line is the best fit 
to the observational compilation of [56]. 




FIGURE 5. Number counts in our model at a variety of wavelengtiis. Tiie top two panels show in two of thie Hubble ACS bands; 
as in previous figures the solid black line represents the fiducial model, and the dashed blue line the low model. Red, blue and green 
data is from the compilation by [34], which includes data from the Hubble Ultra-Deep Field. Additional data from SDSS-DR6 is 
provided by [57]. The lower two plots are infrared bands (Spitzer IRAC and MIPS). Data in the IRAC 3.6 flm band is from [24, 58]; 
the MIPS data is from [29, 59, 27]. 



recent measurements of the fundamental cosmological 
parameters. In Figures 1 and 2 we show the present-day 
extragalactic background light and luminosity density for 
two models of galaxy evolution. Our 'fiducial' model is 
based upon a concordance cosmology with Q.^ = 0.3, 
Q.A = 0.7, Ho = 70.0, and as = 0.9. Our 'low' model 
is based upon the best fit values from WMAP3, which 
include a lower Cg of 0.761. This lower normalization of 
the primordial power spectrum leads to delayed structure 
formalion and lower EBL. The most recent data from 
WMAP5 L67J suggests an intermediate value of ag = 
0.817±0.026, so we believe these two models can serve 
as bracketing the reasonable possibilities arising from 
cosmological uncertainty. The S AMs used here are based 
either upon the merger history of the dark matter halos 
in the simulation or else the extended Press-Schechter 
method, in which for each present-day halo, a merger 
'tree' is constructed, using the method similar to that 
described in [68]. The merger trees track the buildup 
of dark matter mass, with junctions representing halo 
mergers which form liirger, virialized halos. The NFW 
profile [69] is used as the initial halo description with the 
concentration determined using a fitting formula based 
on [70] . The model does not account for scatter or merger 
history in determining concentration. These methods do 
allow tidal disruption and destruction of halos in minor 
mergers to be considered; if a halo is destroyed prior 
to merging, its stars join the diffuse stellar component 
around the central galaxy. 

Gas which cools around the halo is assumed to initially 
fall into a thin disk with exponential profile. The scale ra- 
dius of this disk is determined using conservation of an- 
gular momentum and the concentration and baryon frac- 
tion of the disk. Our model computes the cooling time for 
gas based upon density, metalUcity, and the initial virial 
temperature. Our recipe agrees well with gas infall and 
cooling rates from 3-dimensional hydrodynamic simula- 
tions of cold and hot flows [71, 72, 73 j. 

Feedback from supemovae can heat the cold gas reser- 
voir and drive it from the galaxy. This gas will either be 
deposited in the hot reservoir of the galaxy, or returned to 
the IGM, depending on the wind velocity relative to the 
virial velocity of the halo. Ejected gas can cool and return 
to the galaxy on a timescale roughly equal to the dynam- 
ical time of the galactic halo. The model discriminates 
between cold- and hot- mode accretion of gas based on 
the relative values of the cooling time of the gas vs dy- 
namical time of the halo. 

Star-formation in our model occurs in two regimes, 
quiescent star-formation and merger-driven starbursts. 
The contribution from each of these modes can be seen 
in Figure 3. The integrated mass vs. redshift is plotted in 
Figure 4. The quiescent star formation is treated using 
a recipe based on the Schmidt-Kermicutt law [74, 75] 
with slope 1.4 and normaUzation to the Chabrier IMF. 



Star formation is assumed to cut off below a specific 
surface density, giving the exponential galactic disks a 
corresponding radius within which stars are forming. 
Merger-driven bursts are parametrized by the mass ratio 
of the merging pair, the mass in this case being the total 
mass in the inner part of the halo, taken to be twice 
the characteristic NFW scale radius. The burst efficiency 
parameter determines the fraction of cold gas converted 
to stars in the burst, and no burst occurs for a mass 
ratio of less than 1 to 10. The functional form taken 
for this parameter is from [76]. The SFR during the 
burst is proportional to the available fuel, and therefore 
takes the form of a decaying exponential after an event. 
Details of the functional forms for the burst efficiency 
and timescale can be found in [52]. 

Our implementation of AGN feedback is similar in 
many respects to that of [77], in which AGN operate in 
'bright' and 'radio' modes, with the latter being switched 
on when the accretion rate rises above a critical value. 
All bright mode accretion is triggered by galaxy mergers, 
and the interdependent processes of AGN activity and 
black hole growth are based upon results from a large 
suite of hydrodynamic simulations. Every top-level halo 
in our simulation begins with a seed black hole of 100 
solar masses; we do not find results to be sensitive to 
this amount. Black holes in merging galaxies coalesce 
rapidly and the product grows at the Eddington rate until 
reaching a critical 'blowout' mass, after which accretion 
rate falls as a power law, until it is cut off completely 
upon attaining a final mass determined by the spheroid 
mass and gas fraction. Radiative momentum from the 
accreting black hole is transferred to the galactic wind 
with an assumed coupling efficiency. In radio mode, the 
black hole enters a phase of Bondi-Hoyle accretion [78]. 
The net cooling in this mode is the cooling rate minus 
the heating from Bondi accretion. Heating is ignored if 
the cooling time is shorter than the dynamical time of 
the halo ('cold mode' cooling). 

Predictions for galaxy number counts in the optical 
and IR are shown in Figure 5. The detection of the 
EBL in the far-IR by DIRBE and FIRAS at a level 
comparable to the direct emission from starhght requires 
the inclusion of galaxy populations very different from 
the Milky Way and other nearby galaxies, which emit 
most light in the optical and near-IR. Detectors in the 
mid-IR such as ISOCAM and MIPS find significant IR 
emission from star-forming galaxies. The cause of this 
is primarily tiny particles of dust released in supernovae 
explosions, with typical sizes ranging from nanometers 
to tenths of a micron. The temperature of these grains 
is only weakly affected by the flux of the surrounding 
radiation field, due to the blackbody emission which 
increases as the 4th power of temperature and the falloff 
of emission at long wavelengths, and therefore the peak 
emission ranges from about about 170 microns for a 
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FIGURE 6. The attenuation e^^ of gamma-rays vs. gamma- 
ray energy, for sources at z = 0.05, 0.1, 0.2, and 0.536 (the 
redshift of 3C279). As previously, the fiducial model is solid 
black and the low model is dashed blue. Increasing distance 
causes absorption features to increase in magnitude and appear 
at lower energies. The plateau seen between 1 and 10 TeV 
at low redshift is a product of the mid-IR valley in the EBL 
spectrum. 



Milky Way type spiral to 60 microns for a ULIRG [79]. 
Poly-aromatic hydrocarbons (PAHs) are another class of 
absorbers; tiny particles which emit at a group of specific 
wavelengths in the mid-IR from 3 to 17 microns. The 
absorption and re-emission of light, primarily from the 
UV, by dust and PAH particles in our model is accounted 
for using the STARDUST templates of [63, 64, 80], which 
are embedded in our semi-analytic model, and account 
for emission at wavelengths from the far-UV to sub- 
millimeter These models are based upon the IRAS color- 
luminosity relation, and take the age, star-formation rate 
timescale, and the size of the gaseous disk as inputs. 
The metallicity is tracked both in the stellar spectra and 
surrounding gas, with predictions for radiative transfer 
based on geometry and metal distribution modeled after 
the Milky Way and nearby galaxies. We are working on 
improved dust models that wiU be more self-consistent. 

3. GAMMA-RAY ATTENUATION 

A consequence of the EBL is the attenuation of high- 
energy gamma rays through electron-positron pair pro- 
duction [81]. The threshold for this process occurs when 
the center of mass energy of the photons is equal to the 



rest mass of the electron-positron pair. The cross section 
is maximized at approximately for center of mass ener- 
gies of approximately twice the threshold, and falls as 
inverse energy for E ^ £,/,. Gamma rays above 1 TeV 
are most attenuated by the near- and mid-IR range of the 
EBL. Those in the 300 GeV to 1 TeV regime are sensi- 
tive to light in the near-IR and optical; below 300 GeV 
only UV photons have sufficient energy to cause the pair- 
production interaction. Below 20 GeV there is little ab- 
sorption due to the increasing scarcity of hard UV back- 
ground photons [104]. 

Observations of the attenuation of gamma-ray spectra 
from extragalactic sources provide a measurement of the 
EBL. We predict this attenuation here based on propaga- 
tion of gamma-rays through the evolving EBL given by 
our fiducial and low models. In Figure 6, we show the op- 
tical depth T (attenuation = e^^) for gamma-rays vs. en- 
ergy for a variety of redshifts. In principle, the cosmolog- 
ical history of the EBL could be reconstructed by com- 
paiing observations of high-energy sources at different 
redshifts to their known intrinsic spectra. Unfortunately, 
the emission mechanisms of GeV and TeV sources are 
in general poorly understood. Ground-based gamma-ray 
astronomy, which began largely with the Whipple and 
Hegra collaborations in the 1990s, has opened up the en- 
ergy range above 100 GeV. Today exploration in this 
regime is led by the VERITAS [82], H.E.S.S. [83], and 
MAGIC [84] experiments. The recently-launched Fermi 
(GLAST) and AGILE satellites [85] will provide much 
needed data. The LAT instrument on Fermi has an en- 
ergy range of 20 MeV to 300 GeV, and will efficiently 
scan the sky in this energy range. In addition to detecting 
sources on its own, Fermi will act as a finder for pointed 
observations by ground based telescopes. While in sur- 
vey mode, the LAT will see the entire sky every 3 hours, 
and will be able to report flaring blazars to other detec- 
tors for study. 

Ground-based detectors searching above 100 GeV 
have identified 23 extragalactic sources at the time of 
writing, including 21 BL Lac objects, radio galaxy M87, 
and the flat-spectrum radio quasar 3C279. With the ex- 
ception of M87 these objects are all blazars, accreting 
AGN which generate tightly-beamed relativistic jets that 
lie at a small angle relative to our line of sight. Experi- 
ments at a wide variety of energy ranges have revealed 
that the spectra of blazars consist of a two-peaked dis- 
tribution, the lower peak between the infrared and x- 
rays, the upper in the gamma-rays. This pattern is com- 
monly understood as the result of a relativistic jet of 
charged leptons. The lower peak is attributed to syn- 
chrotron emission, while the gamma-ray peak is the re- 
sult of inverse Compton radiation seeded by either the 
synchrotron photons themselves (the synchrotron-self 
Compton or SSC model) or photons from other sources 
(the external Compton or EC model). It has been ar- 
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FIGURE 7. Here we show limits from gamma-ray experi- 
ments compared to the 4 models presented in Figure 1. These 
limits are based upon the F > 1 .5 criterion discussed in the text. 
The thick green curve is the limit from [86] and the red stars 
are from the observations of the 3C279 by MAGIC [87]. The 
purple points are from functional EBL forms of [88] and the 
point at 60 /im is from [89]. 



gued in the case of at least 2 rapidly varying blazars 
that such variability is better modeled with an external 
radiation source [90]. There is evidence that there exists 
a BL Lac sequence, with less luminous objects having 
both peaks at higher energies than more brilliant sources 
[91, 92, 93], the extreme case being the HBL sources, in 
contrast to the brighter LBL sources. The sequence may 
arise as a result of increased cooling of the accelerated 
particles at higher luminosities, suggesting that sources 
with brighter bolometric luminosities would have lower 
Lorentz factors and peak energies [92, 94]. While they 
account for almost all the detected sources above 100 
GeV, BL Lac objects are themselves only a small subset 
(~20%) of all blazar sources, the other 80 percent being 
flat spectrum radio quasars like 3C279. 

While uncertainties and likely variation in the intrin- 
sic spectrum of blazars make it difficult to link the ob- 
served spectra directly to EBL attenuation, it is possible 
to translate limits on the spectra to EBL constraints. The 
standard assumption in placing limits the EBL from in- 
dividual spectra is that the reconstructed intrinsic spec- 
trum should not have a spectral index harder than 1.5; 
that is, r > 1 .5 if dN /dE ^ E^^ for photon count A^, or 
alternatively dF /dE ^ ii^C"^') for flux F . The standard 
value for a single-zone SSC spectrum is F = (a + l)/2, 
where a is the spectral index of the shock accelerated 



electrons, and is not harder than 1.5 in most models [95]. 
This can be invalidated by assuming a non-standard spec- 
trum for the electrons; a low energy cutoff in the electron 
energy will lead to inverse-Compton accelerated photons 
with an index as low as F = 2/3 [96]. In Figure 7 we 
show a number of constraints on present-day EBL exper- 
iments from gamma-ray observations. The most recent 
limits on the EBL come from observations of blazars at 
more distant redshifts (z>0.1) that have been detected 
by the current generation of ground-based ACTs. Obser- 
vation by H.E.S.S. of two blazars at z=0.165 and 0.186 
were used to set limits on the near-lR EBL based on the 
F > 1 .5 criterion [86]; in this case the maximal Umit was 
our 2001 model [51] multiplied by a factor of 0.45. It 
should be noted that our 2005 EBL prediction [32] in- 
cluded an optical and near-IR flux even lower than this 
older model with the multiplier of 0.45. Another pa- 
per by the H.E.S.S. group set constraints from blazar 
lES 0229H-200 at z=0.1396 [97]. While this blazai- is a 
closer source than the two featured in the 2006 publica- 
tion, the observed spectrum extended above 10 TeV and 
therefore probed the background in the mid-IR. In this 
regime, the effect of optical depth on spectral form is 
minimal due to the approximate falloff in EBL flux. 
The limits derived in this case are well above our 2005 
model for a couple of different spectral slopes consid- 
ered. The most distant source observed at very high en- 
ergies at the time of writing is quasar 3C279 at z=0.536, 
observed by the MAGIC experiment during a flare in 
February 2006 [98]. The spectrum observed was quite 
steep, 4.1 ± 0.7 Stat ± 0.2„ s, and extended from about 80 
to nearly 500 GeV. An analysis of the spectral modifi- 
cation [87] found that there was little room for an EBL 
flux in the optical higher than one consistent with lower 
limits from number counts, approximately equivalent to 
the model of [32] and with our fiducial and low mod- 
els presented here - see Figure 8. An alternative analysis 
of the spectral deconvolution of 3C279 by [99] disputed 
this analysis and argued that the higher EBL of [50, 100] 
could still lead to a steep best-fit spectrum. Another ap- 
proach to the problem is to attempt to constrain the EBL 
by using spectra from several sources simultaneously. 
Dwek and Krennrich [89] considered 12 such permuta- 
tions, and derived an upper limit at 60/im by declaring 
invalid those realizations leading to unphysical intrinsic 
blazar spectra with sharply rising TeV emissions. An ef- 
fort by [101] used observations of 4 blazars by HEGRA 
to limit the EBL in the near-IR and at 60 microns with 
limiting F > 1.0 and 1.5. More recently, this method was 
used in [88], who applied constraints from all observed 
TeV blazars to a large number of possible EBL func- 
tional forms created using a spline interpolation across 
a grid in flux versus wavelength space. The lower bound 
of the union of excluded models formed an envelope rep- 
resenting the highest possible background that does not 



violate any constraints. This was done for 'realistic' and 
'extreme' bounds of F > 1.5 and 2/3 respectively, and 
provided a limit on the EBL from the optical to the far- 
IR. 



4. DISCUSSION 

Figure 8 shows the gamma-ray attenuation edge plot for 
our models. This is a complementary plot to the opti- 
cal depths of Figure 6, and shows the distance to which 
gamma-ray detectors can view before encountering var- 
ious levels of obscuration from photon-photon interac- 
tions. 

New self-consistent hybrid SAMs based on physi- 
cal scaling from numerical simulations and calibrated 
against empirical constraints now enable us to interpret 
the relationship between galaxies, supermassive black 
holes, and active galactic nuclei across cosmic history. 
The latest semi-analytic models are greatly improved 
over the best SAMs that were available four years ago, 
largely because the improving observations of galaxies 
challenged theorists to do better. Although the ACDM 
dark matter backbone for structure formation was al- 
ready well understood, earlier ACDM-based galaxy for- 
mation models suffered from a set of interlinked prob- 
lems - especially overcooling leading to excessive star 
formation in galaxies and clusters, and a resulting failure 
to produce the observed red/blue galaxy color bimodal- 
ity. New SAMs, including the one [52] that our new 
EBL is based on, now account better for the coupUng be- 
tween the supermassive black holes and their host galax- 
ies. "Bright mode" AGN feedback may regulate BH for- 
mation and temporarily quench star formation, but it is 
not a viable long-term mechanism to keep galaxies red 
by preventing them from forming stars, as observations 
require. Low-accretion rate "radio mode" feedback is a 
promising mechanism for counteracting cooling flows 
and preventing star formation in clusters and galaxies 
over long time scales. Such SAMs now accurately pre- 
dict galaxy number counts and luminosity functions in 
all spectral bands and all redshifts except for sub-nmi 
galaxies. They should therefore allow us to predict the 
EBL rather reliably across the wavelengths relevant for 
gamma-ray attenuation, since the far-IR EBL is not rel- 
evant for attenuation of gannma-rays of even the high- 
est observed energies. The predicted range of the EBL 
is consistent with the latest estimates of EBL evolution 
inferred from observations [33]. Thus we have reason to 
hope that this sort of theoretical calculation of the EBL 
and corresponding gamma-ray attenuation is beginning 
to converge. 

There are two main problems. One is that the treatment 
of the crucial near-IR and mid-IR emission of galaxies 
by SAMs is stiU rather primitive. The approach that we 



and others use is based on assuming that most galaxies 
of similar total bolometric luminosity emit radiation with 
similar spectral energy distributions, so that we can use 
"template" SEDs to calculate the IR emission. However, 
Spitzer observations have shown that even nearby this is 
not true, and we have good reasons to expect that it will 
be even less true at higher redshifts. We are therefore 
working on self-consistent models of dust absorption 
and re-radiation, following on [102, 103] which will be 
incorporated into future models of the EBL. 

The other main problem is calculating the ultraviolet 
EBL and the corresponding attenuation of gamma-rays at 
energies below about 100 GeV, which will be observed 
by Fermi (GLAST) and by the new generation of very 
large ground-based atmospheric Cherenkov telescopes. 
Little is known from direct measurements about the EBL 
at energies above the Lyman limit. We do not know the 
relative importance of ionizing radiation from AGN vs. 
stars and the redshift evolution of both. Also, most ioniz- 
ing photons from star-forming galaxies are absorbed by 
local cold gas and dust, with an uncertain fraction es- 
caping to the intergalactic medium. The transmission of 
the IGM is also not fully understood. Predicting optical 
depths to lower energy photons is further complicated 
by the fact that the attenuation edge for an optical depth 
of unity increases to redshifts of several, meaning that 
the evolution of ionizing sources must be understood to 
high redshift. Uncertainty in star-formation rates and ef- 
ficiency, evolution of the quasar luminosity function and 
spectrum, and the possibly changing escape fraction or 
initial mass function all make predicting optical depths 
for gamma rays at high redshift much more difficult than 
studies of local absorption at TeV energies. Preliminary 
calculations are presented by Gilmore et al. [104] in this 
volume, including predictions for gamma-ray attenua- 
tion. 

The Fermi satellite is now performing beautifully, fre- 
quently discovering new, or newly brightened, sources of 
gamma rays. The prospect of many new observations in 
the near future of gamma-rays across a wide range of en- 
ergies from sources at many redshifts means that we can 
expect a great deal of progress as theory and observations 
together clarify the entire spectrum of galactic radiation 
and its sources. 
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